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> ■ ABSTRACT 

m : 

, Aims. The aim is to study the polarization of thermal dust emission to see if the alignment of grain by radiative torques could explain the 
observed relation between the degree of polarization and intensity. Predictions are made for polarimetry observations with the Planck satellite. 
Methods. Our results are based on model clouds derived from MHD simulations of magnetized turbulent flows. The continuum radiative 
transfer problem is solved with Monte Carlo methods in order to estimate the three-dimensional distribution of dust emission and the 
radiation field strength affecting the grain alignment. The minimum grain size aligned by the radiative torques is calculated, and the Rayleigh 
polarization reduction factor, R, is derived for different grain size distributions. We show maps of polarized thermal dust emission that are 
^ predicted by the models. The relation between the intensity and polarization degree is examined in self-gravitating cores. Furthermore, we 
Oh study the effects of wavelength, resolution, and observational noise. 

Results. We are able to reproduce the P/I-relation with the grain alignment by radiative torques. The decrease in intrinsic polarization and total 
^ emission means that sub-mm polarimetry carries only little information about the magnetic fields in dense cores with high visual extinction. 
^ The interpretation of the observations will be further complicated by the unknown magnetic field geometry and the fact that what is observed 
as cores may, in fact, be a superposition of several density enhancements. According to our calculations, Planck will be able to map dust 
polarization reliably when Ay exceeds ~ 2™ at spatial resolution of ~ 15'. 
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1. Introduction 

Magnetic fields play an important role in the dynamics of 
molecular clouds and in the star formation process. In proto- 
stellar cores the topology of the magnetic fields can be pre- 
dicted with models of core formation and evolution. However, 
observations of the magnetic field structure is needed to test 
such models. The best technique for studying magnetic fields 
in molecular clouds is to observe the polarization of the light 
from background stars and the polarization of the thermal dust 
emission, both arising because dust grains are aligned by the 
magnetic field. 

Although the polarization of starlight by aligned interstellar 
dust grains was discovered over half a century ago (Hall . 1949. 
Hiltner n'949> . the processes responsible for the alignment are 
still under study. In a recent review by Lazarian ( I2003t the var- 
ious proposed mechanisms are discussed extensively. 

One of the first processes to be suggested was the paramag- 
netic mechanism (Davis & Greenstein ll95H . which is based on 
direct interaction of rotating grains with the interstellar mag- 
netic field. While this mechanism requires stronger magnetic 



fields than have been uncovered by other means, it may con- 
tribute at very small grain sizes. Purcell ( 1979 1 introduced sev- 
eral processes to make grains very fast "suprathermal" rotators 
and identified one of them, catalytic sites ejecting H2 as the ma- 
jor cause of fast grain rotation. However, it couldn't explain, for 
example, why the observations indicated that small grains are 
less aligned than large ones. 

Fairly recently, the Purcell mechanism was re-evaluated in 
a series of papers by Lazarian and colleagues, pointing out pro- 
cesses such as internal grain wobbling (e.g. Jones & Spitzer 
1967 i Lazarian 1994 Lazarian & Roberge [T997l . grain flipping 
(Lazarian & Draine .l999j and "nuclear relaxation" (Lazarian 
& Draine 1999|l. All these contribute to making the Purcell 
mechanism inefficient. 

The radiative torque mechanism was introduced by 
Dolginov ( 1972.1 and Dolginov & Mytrophanov (1976). Draine 
& Weingartner ( 1996) demonstrated the efficiency of radiative 
torques using numerical simulations, while Abbas et al. ( I2004t 
demonstrated it in a laboratory setup. The radiative torques use 
the interaction of radiation with a grain to spin it up. The pre- 
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dictions of the radiative torque mechanism are roughly consis- 
tent with observations (e.g. Lazarian et al. 1 19971 Hildebrand et 
al. 170001 1. 

In view of this, the radiative torque mechanism is a strong 
candidate as the primary mechanism for grain aUgnment inside 
molecular clouds. Cho & Lazarian ( I2005t showed that even 
deep inside GMCs {Ay < 10) large grains can be aligned by ra- 
diative torques, using a simple, radially symmetric core model. 
Thus, far-infrared and submillimeter polarimetry could reliably 
reflect the structure of magnetic fields deep inside molecular 
clouds. 

Padoan et al. ( 1200 U used numerical MHD simulations 
to model a molecular cloud containing protostellar cores, 
and computed dust continuum polarization for self-gravitating 
cores. They made the simplifying assumption that the temper- 
ature of the dust is uniform, which allowed them to calculate 
the polarization based on gas density and magnetic field vec- 
tors. They calculated the polarization maps for two cases: one 
where the grain alignment mechanism efficiency would be in- 
dependent of the visual extinction and another where it was 
assumed that the mechanism would stop working at a sharp 
cut-off of Ay = 3™, and thus dense cores would not contribute 
to the polarization. The latter case reproduced the observed re- 
sults of decreasing polarization degree with increasing inten- 
sity in the cores, which led them to conclude that the grains 
at high visual extinction appear not to be aligned. This would 
mean that submillimeter polarization maps carry little informa- 
tion on the magnetic fields inside protostellar cores at visual 
extinction larger than Ay «; 3™. 

Like Padoan et al. KOOIV we base our study on three- 
dimensional simulations of magnetized gas flow, but we refine 
the analysis method used to predict the polarized dust emission. 
Radiative transfer calculations are used to estimate the actual 
three-dimensional distribution of dust emission. This naturally 
reduces the weight of high column density regions, because 
there the emission will be reduced by the attenuation of the 
heating radiation field. Similarly, radiative transfer modelling 
allows us to estimate, at every point of the cloud, the radia- 
tive torque efficiency of Cho & Lazarian ('2005V This way we 
can test the hypothesis, that reduced polarization of dense cores 
could be due to the non-alignment of smaller grains. We study 
the general morphology of the resulting dust polarization maps. 
In particular, we look at the relation between total and polar- 
ized dust emission in dense cores. Finally, some predictions are 
made concerning the dust polarization that will be observable 
with the Planck satellite. 

2. Polarization 

2.1. Rayleigh Polarization Reduction Factor 

The Rayleigh polarization reduction factor R is a measure of 
imperfect alignment of the dust grains with respect to the mag- 
netic field (Greenberg 1968; see also Lee & Draine 1985 ) . The 
degree of polarization is decreased, when part of the grains are 
misaligned with respect to the magnetic field. In our case, small 
grains are not aligned by radiative torques while the larger ones 
are. Cho and Lazarian (.2005J give a fitting formula for the min- 



imum aligned grain size in the case of a spherically symmetric 
dark cloud: 

flaig = (log nH)'(Av,iD + 5)/2800 //m, (1) 

where tin is the hydrogen density and Av,id is the visual ex- 
tinction used by Cho and Lazarian (2005 1. They used the ra- 
diation field given in Mathis et al. ( 1983 1, which was given as 
a function of the visual extinction measured along a perpen- 
dicular sightline from the surface of an opaque (in the cen- 
tre Ay ~ 200™), spherical cloud. Thus Ay m measures the 
amount of radiation reaching the grain (see Section 3.2 for 
more). Figure^shows Oaig as a function of Ay iu for four den- 
sities nu- 
cha and Lazarian (12005 > note that their goal is only to ob- 
tain a crude estimate. They make a couple of simplifications 
while deriving the above formula. First of all, they consider the 
anisotropy of the radiation only along the magnetic field, and 
neglect the isotropic radiation. They use the anisotropy factor 
y - 0.7 for their dark cloud, in accordance with Draine and 
Weingartner (I1996> . Our model cloud, as later shown, is much 
different in structure, and has large volumes of optically thin 
medium which are illuminated more or less isotropically. By 
using Eq.[2in those regions, we overestimate the grain align- 
ment efficiency. However, as our main goal is to examine the 
cores within the larger cloud, for which the anisotropy factor 
would be valid, we decide to use Eq. [l] We will return to the 
topic of anisotropy in Section 5.1. The relative attenuation of 
radiation field at different wavelengths, and thus Eq.^ depends 
on the grain size distribution. The study of those effects is be- 
yond the scope of this paper, and Eq.[n should be considered, 
as said before, just an approximation. 
The polarization reduction factor is 

X^"" Ci-ann(fl) da 

where dan is the polarization cross section of a randomly 
aligned grain, n{a) the grain number density, a the grain size, 
flmin the minimum size of the grains, Omax the maximum size, 
and flaig the minimum aligned grain size. 

2.2. Polarized Thermal Dust Emission 

The polarized thermal dust emission is calculated following the 
formalism in Fiege & Pudritz (2000). Self-absorption and scat- 
tering can be neglected at submillimeter wavelengths. 

The Stokes Q and U components are equal to the integrals 

q - J" a f cos 2if/cos^ y ds , (3) 

u - ^ af sin 2\j/ cos^ y ds , (4) 

where or is a coefficient of the particle properties to be defined 
later, / a weighting function related to the intensity of the local 
dust emission, ij/ the angle between the projection of B on the 
plane of the sky and the north, and y the angle between the 
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Fig. 1. Minimum aligned grain size as a function of Ay, id using 
the Weingartner & Draine (2001 1 dust model and Li & Draine 
J2001t grain size distribution. 



local B vector and the plane of the sky. In Padoan et al. 120011 
the weighting function / was the local gas density p, while here 
we consider also the case where the weighting function is the 
true dust emission intensity at the given wavelength, obtained 
from the radiative transfer calculations. 
The polarization angle x is given by 



tan2;i' = -. 
1 

and the degree of polarization P is 



P = 
with 

and 

^2 = 



f ds , 



•/(cos y- 



2 



The coefficient a is defined as 

^pol 



RF- 



(5) 



(6) 



(7) 



(8) 



(9) 



where R is the Rayleigh polarization reduction factor due to 
imperfect grain alignment, F is the polarization factor due to 
the turbulent component of the magnetic field, Cpoi is the grain 
polarization cross section and Cran is the average cross section 
of a randomly aUgned grain. For more details, see Lee & Draine 
(fT985t. 



In our study, F - I because the three-dimensional mag- 
netic field is given by the MHD models and we assume that 
the small-scale structure is resolved in the numerical solution. 
The ratio Cpoi/Cran needs to be fixed, because we want to study 
the effects caused by variations of R. In order to avoid unrea- 
sonably high polarization degree, we choose Cpoi/Cran = 0.15. 
This corresponds to the axial ratio of roughly 1.1 (see Fig. 1 in 
Padoan et al. 2001 ). In our study, we use both uniform 7? = 1 as 
well as R which varies according to the local efliciency of the 
radiative torque, as explained in Section 2.1. 



3. Models 

In this section, we describe the MHD models and the radiative 
transfer calculations that are used in the simulation of the dust 
polarization. 

3.1. MHD models 

In this work, we use the three numerical models - A, B, and C - 
which are described in detail in Juvela, Padoan & Nordlund 



J2001I ). They are based on the results of numerical simula- 
tions of highly supersonic magnetohydrodynamic turbulence, 
run on a 128^ computational mesh with periodic boundary con- 
ditions. The initial density and magnetic field were uniform 
and an external random force was applied to drive the turbu- 
lence at roughly constant rms Mach number of the flow. All 
three models have highly supersonic flows with the rms sonic 
Mach number Ms ~ 10. Models B and C are super- Alfvenic, 
with Ma ~ 10, while model A has equipartition of magnetic 
and turbulent kinetic energy. Ma ~ 1, and thus an order of 
magnitude stronger magnetic field than the other two. Model C 
includes also the eff'ect of self-gravity. Details of the numerical 
code used to solve the MHD equations can be found in Padoan 
& Nordlund (.1999J . 

3.2. Radiative transfer calculations 

In our scheme radiative transfer calculations are needed for two 
purposes. First, the radiative transfer is solved in order to com- 
pute the dust temperature and the resulting emission, based on 
a Monte Carlo code discussed in detail elsewhere (Juvela & 
Padoan l7003l Juvela l2005 i. The total intensity of the dust emis- 
sion is used later to analyze the f //-relations. Once the total in- 
tensity and polarization are computed, radiative transfer effects 
can be ignored, because the dust emission is optically thin at 
the wavelengths considered {A > 100 jum). Second, the radia- 
tive transfer is also calculated in order to estimate the three- 
dimensional distribution of the radiation field inside the model 
clouds, which is needed to compute the minimum grain size 
for alignment. This simulation is done in the V-band (0.55//m) 
and the computed ratios between the local and the background 
intensities, ///bg, are transformed into Ay.eff according to the 
formula 



///bg = exp(-1.0S6 X Av.eff). 



(10) 



The values of Av,eff are then used to estimate the minimum size 
of aligned grains according to Eq. ^ However, this Av,eff is 
not the same Av,id as used by Cho & Lazarian COOSI). Ay.efr 
is a measure of the amount of external radiation reaching a 
grain (see Ea.llO>. In the centre of a homogeneous and spher- 
ically symmetric cloud Ay.eft would be equal to Ay. id (half of 
the Ay through the whole cloud). In inhomogeneous clouds 
the radiation propagates preferentially along low Ay sightlines. 
Therefore, the value Ay gff can be significantly smaller than the 
average Ay that would be calculated as an arithmetic average 
over all directions. As Ay id is calculated along the lowest Ay 
sightline (i.e. perpendicular to the surface), the actual amount 
of radiation reaching the cell is less than in the case where Ay _eff 
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would have the same numerical value. Thus, a given value of 
^v.eff needs to be transformed into a lower value of Ayjo in or- 
der to derive the same radiation field. This is done in our case 
by numerically integrating over all directions to derive the ac- 
tual total radiation reaching the cell as a function of Av,id- This 
allows the matching of the corresponding values of Ay, id and 

Av.eff. 

In the calculations, the dust properties correspond to nor- 
mal Milky Way dust with Ry -3.1 (Weingartner & Draine 
[2001 1 and we use the scattering functions calculated by Draine 
<2003> ^ The cloud models are scaled to a linear size L = 6 pc 
and a mean density n - 640 cm"^. The radiative transfer cal- 
culations depend only on column density, which is < N{H) >= 
1.2 X lO^^cm"^, corresponding to < Ay >~ 6™. However, the 
absolute gas density enters the calculations through the density 
dependent polarization reduction factor/? (see next Section). 

4. Results 

4.1. Polarization maps 

We calculate R in the following manner. First, we transform 
Av.etf to the corresponding Av,id and use Av,id to calculate Oaig 
from Eq. ^ for each cell. We set amin = 0.005 fim. We use 
two power-law distributions, n(a) oc a"^-^, with sharp cut-offs 
at flmax = 0.25 jjm and flmax = 1.0 //m. As a third grain size 
distribution we also use the model of Li and Draine f2001t, 
hereafter L&D, which is not a pure power law and has an ef- 
fective cut-off between the previous a,^^^ limits. Eq.|5]is solved 
numerically to obtain the values of the resulting polarization 
reduction factors for a grid of values of ne and Ay, id- In the 
actual cloud model, the 7?-values are obtained for each cell by 
interpolating from the pre-calculated grid, rather than by solv- 
ing Eq. |2for each cell separately. The polarization reduction 
factor R as a function of AyjD are shown in Fig.|3for the dif- 
ferent grain size distributions and different densities. As the 
maximum grain size increases, the R factor becomes less de- 
pendent on both the density and the value of the extinction, 
Ay,iD- This is the direct consequence of Eq.|2 as a.^\„ reaches 
flmax less rapidly. 

In the calculations of the polarization maps, we use model 
C as our baseline, which is scaled to the size L = 6 pc and 
mean density n = 640 cm"^. This results in an optically thick 
cloud with average visual extinction < Ay >~ 6.3™ through the 
cloud. We also choose 353 GHz as our base frequency. Other 
frequencies and model clouds are examined later 

Figure |3] shows the calculated polarization maps for 
353 GHz. The polarization is calculated using the actual three- 
dimensional distribution of dust emission, that is using the local 
dust emission at the given frequency as the weighting function 
/ (see Eq.|5}- In the left frame the polarization reduction factor 
R varies according to the local values of density and AyjD, as 
shown in the middle frame of Fig. |2] In the right hand frame a 
uniform value of R - 1 was assumed. The difference in the po- 
larization degrees is clearly visible. In particular, with variable 
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Fig. 4. Polarization degree for the whole 128^ pixel maps, 
viewed from the y-direction. The uppermost frame corresponds 
to model A where the variable R factors are calculated using Li 
& Draine (2001 1 grain size distribution. The other frames are 
for model C with uniform R - 1 (middle frame) and variable R 
(lower frame), respectively. 



R the polarization degree decreases in regions of high intensity, 
while no such trend exists in the right hand frame. 

In Fig. 13 the polarization degree and intensities are shown 
for the whole maps. The shapes of the plots are quite similar for 
the two lowermost frames, which are for model C with constant 
and variable R factors. The wide range of polarization degrees 
at a any given intensity is caused by the geometry of the mag- 
netic field and not by any numerical noise. This can be clearly 
seen when comparing the results with the uppermost frame, 
which is for model A. Model A has a strong, structured mag- 
netic field in the z-direction, with only small deviations. Thus, 
the depolarization effects of tangled magnetic fields in the line 
of sight are small, leading to a much tighter relation between P 
and /. 

On closer look the two lowermost frames in Fig. |3 show 
some clear differences. First, there is a reduction of the max- 
imum polarization degree from 13% to 9%. This is easily ex- 
plained if one looks at the middle frame of Fig. Hand recalls 
that we are considering a quite dense cloud, leading to R < I 
even in the low-density parts of the cloud. The low-intensity 
points to the extreme left are unaffected, but otherwise the 
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Fig. 2. Polarization reduction factor 7? as a function of Av.id using Weingartner & Draine (200V\ dust model. The curves corre- 
spond to different densities. The frames correspond to different values of the maximum grain size, Omax- In the first and the third 
frame the grain size distributions are pure power laws, n{a) ~ a"^^, while in the middle frame the size distribution is taken from 
Li & Draine (12001^. 
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Fig. 3. Simulated polarization maps at 353 GHz with variable R (left frame) and uniform R - \ (right frame). Polarization vectors 
are drawn for every third pixel, and the scaling of the vector lengths is the same in both frames. In the left frame the maximum 
polarization degree is 9% while in the right frame it is 13%. The background image shows the total intensity at this frequency. 



polarization degree is clearly reduced, even for low intensity 
sightlines. 



Another noticeable feature is the change in the upper enve- 
lope of the polarization degree. This can be seen by comparing 
the lower frames. With variable R, on sightlines with higher 
intensity (and higher average volume density) the maximum 
polarization degree becomes gradually smaller In model C the 
effect is clear, but still small compared with the overall scat- 
ter in the polarization degrees. The overall slope is seen more 
clearly in in the case of model A (the uppermost frame), be- 
cause of the smaller scatter. The maximum polarization drops 
from about 12% to 6% in the interval 0.06 < ///max < 1-0. 
These effects will will be discussed further in the next subsec- 
tion. 



4.2. Self-gravitating cores 

Figures |5] and |6l show zoomed 20 pixel by 20 pixel maps of 
selected cores of model C clumps. Maps are shown for both 
R = I and variable R. The model C included the effect of self- 
gravity. However, what is interpreted as a core in observations 
need not be a clear-cut, compact gravity bound cores in reality. 
The same applies to these objects that were selected based on 
observations in the direction y. These are knots in the filamen- 
tary cloud structure, and in certain directions the structures can 
be very elongated. 

The general trend towards smaller polarization degree in 
the cores is obvious. The relation between observed total in- 
tensity and polarization degree is, however, not always clear 
Density enhancements close to cloud edges provide one excep- 
tion. Because they are subjected to unattenuated external radi- 
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ation field, they can exhibit large emission although they have 
relatively low density compared with actual cores. Thus, the 
factor R is larger than in the cores and the polarization very sim- 
ilar to the case of uniform R. One example of this can be seen in 
Core 3, when Figs.|5]and |6lare compared. This is particularly 
striking in the direction x, where the actual core is situated in 
an area of lower brightness slightly above the map centre, while 
a bright filament on the cloud surface crosses the lower part of 
the map. Another thing to note is that the reduction of the po- 
larization in the dense parts of the cloud can, in special cases, 
actually increase the total polarization degree. The core 2 pro- 
vides an example of this, when observed from the x-direction. 
Although the polarization decreases in the bright core itself, the 
polarization increases immediately to the left of it. This is, of 
course, due to the magnetic field geometry. If the magnetic field 
is predominantly of opposite direction in high and low density 
regions on the line of sight, they depolarize the observed inten- 
sity. When the polarization reduction factor reduces the polar- 
ization in the dense region, the depolarization effect is reduced 
and the observed net polarization is correspondingly increased. 
This is another indication of the fact that it is difficult to inter- 
pret observations of polarization if not in a statistical sense. 

Fig. shows the plotted f //-relations for the cores seen 
in Fig. |6l In the direction x plot of Core 3, one can clearly 
identify the bright, high polarization filament on the edge of 
our cloud. Other bright but less dense (thus high polarization) 
regions help to mask the slope. The observed polarization is 
a sum of many effects, depolarization being the most obvious 
one. For example, the large-scale structure of the cloud deter- 
mines which parts of the cloud receive the most radiation from 
the external radiation field. This affects the efficiency of the 
radiative torque and of grain alignment (taken into account in 
'^v.iD in Eq. [0 and also heats up the dust grains, which in- 
creases their thermal emission. In the plots, the scaling of the 
.X-axis depends on the maximum intensity which, therefore, af- 
fects the appearance of the figures. One example of this is Core 
2 seen in the direction y, where one can easily see the bright 
upper right corner of Fig. |6l represented as the cluster of high 
polarization, high intensity points in Fig. If this area were 
omitted, the drop of polarization degree with intensity would 
become much more dramatic. 

In Fig. Owe look at some of these effects. In the upper- 
most frame we show for the Core 2 the polarization degree in 
the direction z as a function of column density A^. The clump 
of low polarization points, identified with the filament in the 
lower portion of the image in Fig.|6l moves from moderate in- 
tensity to low column density. Thus these points are not part 
of any high-density core, and the low polarization is likely to 
be caused by the depolarization resulting from the magnetic 
field structure. If those points were omitted, we would end up 
with a quite clean P//-relation for the core itself. Of course, the 
problem with this kind of approach is that, for an observer, the 
true column density is not a priori known. However, if accu- 
rate column density estimates were made, this kind of point-to- 
point investigation would help in distinguishing between true 
cores and smaller density enhancements on the cloud surface. 
If the magnetic field is preferentially parallel to the filament, 
the result, high intensity and low polarization, mimics the ap- 
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Fig. 8. Selected cases for comparison with Figs. |5l and |6l The 
upper frame shows the polarization degree for Core 2 as a func- 
tion of true column density, when viewed from the direction z- 
The lower frame shows Core 1 viewed from the direction x, 
with variable R (dots) and uniform R - 1 (crosses). The polar- 
ization degree is plotted as a function of the total intensity at 
the same frequency of 353 GHz. 



pearance of a dense core with small value of R. There are some 
example of the opposite effect in the cores, where the fainter 
core points move in the plots to the right when we shift from 
the intensity to the column density axis. Long filaments that are 
aligned along the line-of-sight would produce a relatively high 
surface brightness compared with the actual column density. 

In the lowermost frame of Fig.|S]the difference between the 
uniform R = I case and the variable R depending on density 
and Av,iD is obvious. With R = I the maximum polarization 
increases, rather than decreases, to ~IQ% with increasining in- 
tensity. When R is variable there is a downward trend, ending 
at around 3% at the brightest regions. This is explained by the 
fact that intensity and density are usually connected. The high- 
intensity points correspond to the high column density. In the 
high density regions we also expect to find high Ay, id- Thus, 
looking at the middle frame of Fig 121 we see that the polar- 
ization reduction factor R decreases quite rapidly as we go to 
denser and more opaque regions. Therefore, those regions af- 
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Fig. 5. Selected cores of model C observed at 353 GHz. The background grayscale map shows the intensity, and the vectors the 
direction and degree of polarization for the case of 7? = 1 . 



feet the polarization less, whereas in the R - I case, it is pre- 
cisely the dense parts that contribute the most. 

4.3. Effects of Different Cloud and Dust Models 

In this section, we examine differences between diiferent cloud 
models as well as discuss the possible effects of grain growth. 
We will also briefly consider the effect of the choice of wave- 
length in our polarization calculation. 

Figure. |5] shows some histograms for out three model 
clouds. Particularly noticeable is the sharp contrast when the 
polarization degree and angle of model A are compared with 
the other models B and C. Similarly, in model A there is a very 
strong dependence on the direction of observations. These re- 
sults from the fact that model A has a very strong magnetic field 
that is aligned along the z-direction. There is Uttle depolariza- 



tion caused by magnetic field tangling, and, consequently, the 
polarization degree is high, and the angle of the polarization 
vectors is almost uniform. When model A is observed from the 
z-direction, the angle y is close to perpendicular to the plane 
of the sky and hence cos^ y is small (see Section 2.2). Thus, 
the average polarized signal is small, and small variation in 
field geometry can lead to strong depolarization. However, be- 
cause there is an obvious preferred angle even in z-direction, 
we can say that the evolved magnetic field is not aligned per- 
fectly along the z-axis, but has a small deviation from it. 

Models B and C are both super- Alfvenic and, consequently, 
the field geometry is much more complicated. Comparing these 
two models, we notice that the polarization degree is slightly 
smaller in model C when observed from any direction. On the 
other hand, the distributions of polarization degrees are quite 
similar. While model C has slightly more very dense cores. 
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Fig. 6. Selected cores of model C at 353 GHz. The figure is identical with Fig.|5l but assumes variable polarization reduction 
factor, R. 



probably due to its inclusion of self-gravity, the number of 
dense cells is so small that it seems implausible that R would be 
the cause of this difference, especially since model B shows a 
similarly large fraction of very dense cells. However, model B 
seems to have a slightly less tangled magnetic field and, there- 
fore, less depolarization. 

In Fig. no|we show similar histograms for model C when 
the cloud is observed at resolutions corresponding to 1,4, or 
16 pixels. The four-pixel resolution seems almost identical to 
the original one. The 16-pixel resolution shows clear results of 
the averaging, by narrowing down the variance of the distri- 
butions. Of course, this is to be expected as soon as the beam 
size exceeds the size of typical structures. If we are interested 
in the details, such as the cores, the best possible resolution is, 
of course, needed. On the other hand, for general morphologi- 



cal and statistical studies the results are not very sensitive to a 
small reduction in resolution. 

Figures ITTI and [21 are similar to Figs. |6l and showing 
Core 2 in y-direction with the three different grain size distribu- 
tions that were introduced at the beginning of Section im The 
increase of the polarization vectors and even a slight change of 
direction is clearly visible, as we move toward larger amax- This 
can be easily explained by looking at Fig|2 The smaller the pa- 
rameter flmax, the faster R decreases as a function of density and 
Ay. ID- This is because in Eq.|2]flaig quickly reaches the smallest 
flmax, while the largest Omax takes a larger range of Ay, id values 
to build up to, hence making the transition more smooth. 

In Fig.[21we see the drop in the maximum P, similar to the 
global plots. One can also see that while the two frames on the 
left resemble one another in shape, the rightmost frame is sig- 
nificantly different. This hints at the fact that at Omax - 1 jum. 
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Fig. 7. The relation between polarization degree and total intensity for selected cores of the model C. The polarization degree 
has been calculated using the WD2001 dust model and variable polarization reduction factor, R. 



we are approaching the case of R - 1. This indicates that the 
number of very dense cells is relatively small, as can be seen 
from Fig. |9l and that Av.id does not reach as high values as 
Av, the visual extinction through the whole cloud. In an in- 
homogeneous cloud there are always low density regions and 
'holes' that allow external radiation field to penetrate deep into 
the cloud. Therefore, Av,id values can be much smaller than 
the average Ay, which, furthermore, is calculated through the 
whole cloud. In model C, Av,id has a maximum value of 4.3™ 
and a mean value of 0.03™. We can also see the slope of the de- 
crease of maximum P slowly leveling off as we move to higher 

values of Omax- 

In Fig.^lwe have calculate the polarization for Core 2 un- 
der different assumptions. In the first frame we use the Padoan 
et al. 1200 II approach, and assume that the dust emission is di- 
rectly proportional to the density. In practice, this means that 
the weighting function / in Eq.|3]is simply the local density. 
As a further consequence, the selected frequency has no ef- 
fect on the polarization and affects only the intensity. The other 
two plots are calculated by taking into account the actual three- 
dimensional distribution of dust emission, i.e., the weighting 
function / is the local emission at the frequency in question. 
On the other hand, in these plots the polarization reduction fac- 
tor R was assumed to be constant, R - I. There is a clear, 
systematic change whereby the dynamic scale of the intensity 



axis becomes more narrow. The overall shape of the figures 
remains rather similar, and there is only a small reduction in 
the expected polarization degree. The figures do, however, indi- 
cate that a change in the wavelength affects the relative weights 
that are given to different regions. When one moves to shorter 
wavelengths, the correlation between column density and ob- 
served intensity decreases. Optically thick regions contain cold 
dust and, therefore, radiate mainly at longer wavelengths. If 
one goes down to far-infrared (~ 200yum or below), the emis- 
sion comes preferentially from optically thin regions and from 
the surface of optically thick cores and filaments. Therefore, 
shorter wavelengths work in a similar way as the variable po- 
larization reduction factor, and the observed polarization map 
contains less and less information about magnetic fields inside 
the cores. Thus, the dense cores become indistinguishable from 
their surroundings. 

4.4. Predictions for the Pianck satellite 

Planck^ is a CMB mission of the European Space Agency 
(ESA) that is planned to be launched in 2007. The satellite wiU 
provide sensitive all sky maps in a wide frequency range be- 
tween 30 and 857 GHz. In the shortest wavelength channels 
the Galactic dust emission will be the dominant emission com- 

- http;//www.rssd.esa.int/Planck/ 
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Fig. 10. Histograms of column density, polarization degree, and polarization angle computed for model C viewed from y- 
direction. The three curves correspond to different resolution of observation, corresponding to FWHM of either 1, 4, or 16 pixels. 



ponent. Planck has polarization capabilities in seven channels 
(see Table^. In the following we examine the observable dust 
polarization at two frequencies, 143 GHz and 353 GHz. We use 
monochromatic simulations and no convolution is done over 
the band profiles. However, the frequencies studied are rela- 
tively far in the Rayleigh-Jeans regime, so that this simplifica- 
tion is not expected to have a large effect, apart from a possible 
small error in the estimated absolute values of the intensities. 
Our previous tests showed already that morphological differ- 



ences will arise only in comparison with much shorter wave- 
lengths (see Fig.ll3>. 

For this test we use model C. The model was previously 
scaled to represent a cloud with a size of 6pc and mean den- 
sity 640 cm"^, resulting in an optically thick cloud with aver- 
age visual extinction < Ay >~ 6.3™. This was done in order 
to bring out the effects associated with high densities and high 
values of the effective extinction (e.g., in Eq.^. Here we want 
to study whether polarized signal might still be detectable in 
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Fig. 13. Polarization degree of Core 2 viewed from the y-direction. The size distribution of the grains is taken from Li & Draine 
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Table 1. The performance goal values for the resolution and sensitivity of the Planck satellite. The ex- 
pected sensitivities after 14 months of observations (two full sky surveys) are given as thermodynamic tempera- 
ture relative to the CMB temperature of 2.73 K. A more complete table is available through the ESA website 
|http : //www . rssd . esa . int/?pro j ect=PLANCK?(page=per£_top 



Centre frequency 
(GHz) 


Bandwidth 
(GHz) 


Beamsize 
FWHM (arcmin) 


Sensitivity, / 
[AT/TxlO-*] 


Sensitivity, Q and U 
[AT/TxlO-"] 


Sensitivity, / 
[10' Jy/sr] 


Sensitivity, Q and U 
[10^ Jy/sr] 
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regions of low column density. Therefore, we scale the column 
densities down by a factor of ten, resulting in a model with av- 
erage extinction < Ay >~ 0.6™. Simulations are done using a 
variable polarization reduction factor R that corresponds to the 
Weingartner & Draine J2001> dust model. Noise is added to the 
total intensity and the Q- and U-components according to the 
values given in Table [2 

Figure shows the computed polarization map for 
143 GHz. In the figure we have also plotted contours cor- 
responding to Ay values of 0.5 and 2.5 magnitudes. Below 
Ay ~ 0.5 the polarized signal is dominated by noise, which 
randomizes the angles and increases the average polarization 
degree. Only above Ay ~2.5™ is the signal comparable to the 
noise. On the other hand, for total intensity the signal to noise 
ratio is high and the noise has no visible effect on that pa- 
rameter The average intensity over this map was 1.8x10"* Jy/sr 
while the 1 - cr noise, 2.3x10^ Jy/sr, is smaller about one order 
of magnitude. 

In Fig.^Ja similar plot is shown for 353 GHz. In terms 
of absolute intensity, the sensitivity is lower at this frequency. 
However, the rise in the dust spectral energy distribution more 
than compensates this, and the signal-to-noise ratio is signifi- 
cantly higher. The polarization vectors follow, with some scat- 
ter, the magnetic fields even along filaments with Ay ~ 1™ or 
more. At 353 GHz the spatial resolution of Planck is 5 arcmin 
so that the pixel size is somewhat smaller than in Fig. 1141 

The sensitivity can be increased by reducing the spatial res- 
olution. In Fig. ^] the pixel size of the previous maps is set 
equal to 5 arcmin, and the maps have been convolved down to 
a resolution of FWHM - 15 arc minutes. The signal-to-noise 
ratio is, of course, correspondingly larger by a factor of ~ 3. At 
353 GHz a comparison with Fig.ll5h shows that the orientation 
of the polarization vectors is generally correct above Ay ~ 1", 
and is no longer completely random even in the regions with 
the lowest column density. At 143 GHz the scatter is larger, and 
thus even in the regions with Ay ~ 2.5™ the true polarization 
field cannot be reliably measured. 

The errors of polarization degree P (solid line, left hand 
scale) and polarization angle (dashed line, right hand scale) of 
the convolved 353 GHz map are shown in Fig. El At Ay ~ 1™ 
the relative error of P is ~40% and x ~30 degrees. At Ay ~ 



2™ the errors are already relatively low, ~25% for P and ~15 
degrees for;^'. 

The 143 GHz band is already close to the peak of the CMB 
spectrum but, even at these column densities, the total signal is 
dominated by the dust emission. In the previous 143 GHz map 
the minimum column density is l.lxlO^'^cm"-', corresponding 
to less than 0.1 magnitudes of visual extinction and a surface 
brightness of ~ 3 x 10^ Jy/sr At this frequency, the intensity 
of the CMB signal, 3.8 x 10^ Jy/sr, is still lower by almost one 
order of magnitude. Because the polarization degree is of same 
order for both the dust emission and the CMB, also the polar- 
ized signal will be dominated by the dust emission. Conversely, 
the CMB signal introduces errors to the estimated dust polar- 
ization signal only at the 10% level. If one wishes to study the 
dust emission at any fainter level, the CMB should be taken into 
account as a possible error source. Of course, at higher frequen- 
cies the contribution of the CMB signal becomes negligible, 
and the combination of several wavelengths makes dust polar- 
ization studies possible even at fainter levels, if the resolution 
is correspondingly decreased. Conversely, in Planck CMB ob- 
servations, the dust contamination is minimized by using mea- 
surements at high galactic latitudes and lower frequency bands, 
70 GHz and 100 GHz. 

5. Discussion 

5.1. Current polarization model 

We have shown dust polarization maps computed for model 
clouds that are based on three-dimensional MHD simulations. 
The main improvement with respect to the previous study by 
Padoan et al. J2001> is that we have calculated the true distri- 
bution of dust emission with the aid of three dimensional ra- 
diative transfer calculations. This affects the way dilferent vol- 
ume elements contribute to the observed signal. Furthermore, 
instead of resorting to an ad hoc model, the polarization reduc- 
tion factor has been calculated according to the model of Cho 
& Lazarian J2005> . In this model the grain alignment depends 
on radiative torques and, consequently, radiative transfer calcu- 
lations have also been used to estimate the strength of the ex- 
ternal radiation field at each position of the model clouds. This 
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Fig. 14. Simulated polarization map at 143 GHz without noise 
(upper frame) and with noise corresponding to observations 
with the Planck satellite (lower frame). The cloud is model C, 
scaled down to < Ay >=0.6'". One pixel corresponds to 7 arc 
minutes, the Planck resolution at this frequency. Polarization 
vectors have been drawn for every third pixel, and the scaling 
of the vectors is the same in both frames. In the upper frame 
the maximum polarization degree is 9%. In the bottom frame 
noise dominates and, for clarity, long vectors are cut at 25% 
and replaced by a dot. The background image shows the to- 
tal intensity, and the contours denote Ay values of 0.5 and 2.5 
magnitudes. 

has been accomplished through the determination of values of 

Av.eff- 

Observations show a decreasing polarization degree toward 
dense cores. Gon§alves et al. ( I2005> proposed a geometrical ef- 
fect, the toroidal magnetic fieldlines bending to resist the grav- 
itational pull, to partly explain this drop. Padoan et al. (.2001 ) 
could explain the observed reduction by a qualitative model 
where the dust alignment was restricted to regions with Ay be- 
low 3". This extinction was calculated in similar fashion to 
our Ay.eif . With our models we find a very similar behaviour in 
the plots of polarization degree vs. intensity (e.g.. Fig. |6}. At 
low column densities the polarization degree varies between 
~4 and 10% while, at the centre of the cores, the degree drops 




Fig. 15. Simulated polarization map at 353 GHz without noise 
(upper frame) and with noise corresponding to observations 
with the Planck satellite (bottom frame). Without noise the 
maximum polarization degree is 9%. In the upper frame map, 
long vectors are cut at 25% and replaced by a dot. One 
pixel corresponds to the Planck resolution of 5 arc minutes. 
Polarization vectors are drawn for every third pixel along both 
axes. 

to a few percent or, in some cases, even below one per cent. 
In other words, our physical model is capable of explaining 
the general observed behaviour (e.g., Ward-Thompson et al. 
12000, Henning et al. 1200 It . At the same time, the individual 
fields show considerable variations, and we find also cases with 
no apparent decrease in the polarization degree (e.g.. Core 3 
and direction z) or even increased polarization around intensity 
maximum (Core 3 and the filament seen toward direction x). 
This indicates that observations of individual cores may show 
considerable variation, depending on the actual configuration 
of the magnetic field. One apparent problem is that from ob- 
servations it may be difficult to discern real and possibly grav- 
itationally bound cores from projection effects. In particular, a 
core with large polarization reduction factor may appear very 
similar to a filament with major axis and field lines parallel to 
the line of sight. 

The individual contributions of local density and Ay.iD to 
the reduction of the polarization degree also merit a mention. In 
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Fig. 16. The effect of convolution on the polarization maps 
shown in Fig. [21 and [21 One pixel of the original simulation 
is assumed to correspond to 5 arc minutes, and the maps are 
convolved down to a resolution of 15 arc minutes. 




Fig. 17. The errors of polarization degree P (solid line, left 
hand scale) and polarization angle x (dashed line, right hand 
scale) of the convolved 353 GHz map in Fig. 1161 



Model C the local density varied between ~ 1 and ~ 10^ cm"^ 
while Ay, ID had values from 0.005 to 4.3. Just by looking at 
Eq.[nit is easy to see that in global scale, the variation in lo- 
cal density dominates over Ay, id- However, in the cores, where 
the densities are around the same order of magnitude, the con- 
tribution of Ay ID comes to its own by emphasizing the differ- 



ence between the illuminated regions on the outer parts of the 
core and the shadowed, even denser interior In cores, both lo- 
cal density and AyjD are of similar importance. 

Related to this issue is the isotropic and anisotropic com- 
ponents of the radiation field, which play major part in the effi- 
ciency of grain alignment by radiative torques. In this work, we 
neglected detailed simulation of the anisotropy of the radiation 
field needed to derive the actual values of the anisotropy fac- 
tor, y, and instead used the fitting formula (Eq. by Cho and 
Lazarian ( 2005 ) where y - 0.7. Draine and Weingartner ( 1996 1 
give 7 = 0.1 for diffuse clouds, due to the anisotropy of the star 
distribution in the Galaxy. We can qualitatively estimate the ef- 
fect that the change of y would have as it would be similar to 
the increase of local density by the corresponding inverse fac- 
tor. For example, the difference between 10^ and 10^ cm~^ at 
low Ay ID (see Fig.|3 is a drop from 0.95 to 0.8 in R. Thus, we 
would expect a slightly lower polarization in those regions, but 
not by much. If anything, it would help to reduce the 'contami- 
nation' of the diffuse background and make the cores stand out 
more. Also, because we would expect to see a more isotropic 
field at the centre of the core than at its exterior, the inclusion of 
y would make the core centres drop in polarization even more 
sharply. In other words, our model is likely to overestimate the 
polarization at the core centres, rather than the opposite. 

In the modelling we assumed dust properties to remain con- 
stant throughout the model clouds, while only the alignment of 
the grains was dependent on the local density and the visual 
extinction. However, there is ample evidence that optical dust 
properties do change already above a few magnitudes of Ay. 
If this were caused primarily by a change in grain size distri- 
bution, this could counteract the reduction obtained via the R 
factor. As an example, let us consider an optically thick core 
with density n ~ 10^ cm""* and AyjD ~ 6". According to Fig.|2l 
the /^-factor could increase by a factor of several if the L&D 
grain distribution were replaced with a distribution with much 
larger flmax- On the other hand. Fig. 1121 shows that, in the ac- 
tual model, the effect of an increasing Omax was rather limited. 
In that particular case, when L&D dust model is replaced with 
MRN type distribution with Omax ~ 1/um, the polarization de- 
gree increases in the core centre only from ~2% to ~3%. It 
is very unlikely that in such a core a„iax could grow close to 
one micron. Another explanation for the relative insensitivity 
to flamax is that the observed intensity is weighted toward low 
density material. There are two reasons for this. First, the polar- 
ization reduction factor decreases rapidly, whatever the value of 
flmax- Second, in a dense core the dust is colder and, therefore, 
the emission is lower for a given column density. In a core with 
Ay.iD ~ 7™ to the centre, the dust temperature may already be 
reduced to ~ 14K or below (Stamatellos & Whitworth 2003}. 
This can be compared with ~18K in diffuse material. In the 
Rayleigh-Jeans tail the effect is small, and even at 353 GHz 
only a factor of 1.5. However, the effect becomes important 
when one approaches the peak of dust emission. At lOOyum the 
intensity produced by an 18 K grain is akeady ten times as large 
as the intensity produced by a 14 K grain. From Fig.|2lone can 
see that at low densities the slope of the R vs. Ay id relation is 
very similar for all values of Omax, and the percentile difference 
in the absolute value of R are similarly small. 
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The above picture may change if dust grains undergo more 
profound changes. The sub-mm emissivity of dust has been 
observed to increase significantly already in clouds with vi- 
sual extinction only a few magnitudes (Bernard et al. 1999 
Stepnik et al. 2003 1. This has been interpreted as evidence of 
not only grain growth but the appearance of dust aggregates 
and fluffy grains of low packing density. At long wavelengths 
the high emissivity might again increase the contribution that 
the densest regions have on the total intensity observed towards 
cores. On the other hand, the formation of dust aggregates can 
also affect the intrinsic polarization efficiency of the grains and 
change the efficiency of the alignment mechanisms (Lazarian 
& Efroimsky 1999 Gupta et al. .2006} . Such effects are not 
considered in this paper 

5.2. Comparison with otiier simulations 

On larger scales, Prunet et al. (I1998> presented simulations of 
polarized dust emission at high Galactic latitudes. The calcula- 
tions were based on three-dimensional cloud structure that was 
estimated based on HI data. The HI line velocity provided in- 
formation on the line-of-sight density distribution. The simula- 
tions assumed constant dust properties with intrinsic polariza- 
tion efficiency ~30%. Furthermore, three field geometries were 
considered. A uniform field direction resulted in high polariza- 
tion degree with very narrow distribution, and was dismissed as 
being incompatible with observations. In our simulations this 
corresponds closely to model A, where the field lines have re- 
mained almost parallel because of the strength of the magnetic 
field. Of course, strong polarization also requires that the cloud 
is not observed parallel to the field lines. In our model A, when 
the sightlines were perpendicular to the field lines, the polar- 
ization degree peaked around 10%. For a random orientation 
the mean polarization would be lower by half. In the other two 
cases studied by Prunet et al. the field was assumed to be either 
parallel with the major axis of the cloud structures, or to have a 
random orientation perpendicular to it. Both resulted in a wider 
distribution of polarization degrees that, in the case of R = 1, 
peaked around 2-3%. These are again similar to our models B 
and C, where, irrespective of the viewing angle, we have an 
approximately log-normal distribution of polarization degrees 
that peaks around the value of ~3%. Therefore, although the 
results were obtained by entirely different methods, the basic 
properties of the dust polarization seem very similar. This is 
reassuring for attempts to build statistically meaningful all-sky 
templates for the polarized dust emission that will be observed, 
e.g., by the Planck satellite. 

Another important statistic is the power spectrum of the po- 
larized signal. Figure^]shows the power spectra computed for 
the model C. Spectra are shown separately for the I, Q, and U 
components, as well as the polarized intensity Ipoi, from di- 
rection y. These can be compared with the power spectra pre- 
sented by Prunet et al. tl998J . The slope of intensity, in inter- 
val 2 < k < 20, is about k~^, in agreement with Prunet et al. 
J1998> . The Stokes parameters Q and U have a slightly more 
shallow slope, but not as flat as by Prunet et al. Cl998j . It can 
also be seen that in small scales, the power in the polarized 




Fig. 18. The power spectra of the polarized intensity Ipoi and 
the Stokes parameters Itot, Q and U are plotted for model 
C from direction y. Column density follows the intensity 
power spectrum closely, and is thus not shown. While model 
A showed a distinct shift in the level of the Q and U power 
spectrum in direction z, it - and model B - were not different 
enough in shape to justify another figure. 



components drops slower than the total intensity, making the 
polarized power an appreciable fraction of the total intensity 
power The reason for this would appear to be the magnetic 
field morphology causing abrupt, pixel to pixel variations or 
larger structures, to about 10-pixel level, while the background 
intensity stays about the same. 

5.3. Comparison with observations 

Direct observations of polarized dust emission are still very 
scarce for all but the very densest regions. The first detection 
of diffuse polarized dust emission was reported by Benoit et 
al. (I2004> . The observations were done with the balloon-borne 
Archeops experiment, with instruments that closely resemble 
those of the Planck satellite. Polarized signal could be detected 
at 353 GHz with the resolution of 13 arcmin, the average po- 
larization degree being of the order of 4-5%. However, a large 
variation was seen between different regions. In the Cygnus 
region the polarization was low, as can be expected if the mag- 
netic field is preferentially parallel to the spiral arms. On the 
other hand, in some regions of the Galactic plane the polar- 
ization reached values of ~15%. In our models, polarization 
degrees this high could be reached only in model A and only 
when the polarization reduction degree R was equal to one (see 
Fig. |4}. However, our maximum polarization degree is some- 
what arbitrarily set by our choice of a, or to be more exact, 
our choice of Cpoi/Cran = 0.15. This ratio scales our maximum 
polarization degree to about 15.3%, which explains why we 
do not reach P - 15% in all our models, since that would 
demand a total alignment and no depolarization on the line 
of sight. A slightly larger, yet still a moderate deviation from 
sphericity, can double our polarization degrees: at axial ratio 
of 1.2, Cpoi/Cran - 0.3. Thus One needs to be very careful not 
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to overstate the case. Certainly a high polarization degree hints 
to an ordered magnetic field and high polarization efficiency 
throughout the cloud, yet even small changes in our assump- 
tions about the grain shapes or lowering the density of our very 
dense model clouds would change the maximum polarization 
degree so that the other two models could duplicate even the 
high polarization degree of 15%. 

Ponthieu et al. ( 12005 > presented the measurements of tem- 
perature and angular power spectra of the diffuse emission of 
Galactic dust at 353 GHz as seen by Archeops on 20% of the 
sky. They concluded that if the fraction of the sky observed by 
Archeops is representative, dust polarized radiation will be a 
major foreground for determining the polarization power spec- 
tra of the CMB at high frequencies above 100 GHz. However, 
their power spectra cannot be directly compared with our re- 
sults, because they cover much larger area, where the turbu- 
lence is no longer the dominant factor affecting the structure of 
the density and the magnetic fields. 

5.4. Future 

In Fig. ^] it is shown that when the spatial resolution of the 
Planck satellite is reduced to 15 arcmin, it will be able to pro- 
vide a relatively clean map of the polarization of a diffuse cloud 
at 353 GHz. The general morphology of the magnetic field is 
more or less correct above Ay ~ 0.5", even if the polarization 
degree P itself has a relative error of ~60%, as seen in Fig. 1171 
The situation improves with increasing Ay, until at Ay ~ 2", 
the errors have dropped to the relatively low values of ~25% for 
P and ~15 degrees for the polarization angle Thus, we con- 
clude that Planck will be able to trace the general morphology 
of the magnetic fields at Ay ~ 1™ and to map dust polarization 
reliably when Ay exceeds ~ 2™. 

In this study we have used simulations with relatively low 
spatial resolution. The results remain reliable at large scales, 
and the discussion of core properties is appropriate for gen- 
eral density enhancements and the initial phases in the forma- 
tion of gravitationally bound cores. In a rotating and collaps- 
ing core the magnetic field lines become wound so that both 
the field strength and the depolarization caused by the field ge- 
ometry increase. If the system were not perfectly symmetric 
one might, at high resolution, detect sharp variation in both the 
strength and direction of the polarization vectors. However, in 
order to be able to follow the development of magnetic fields 
inside such collapsing objects, simulations with considerably 
higher spatial resolution are needed. Spatially variable y and 
the inclusion of other grain alignment mechanisms would al- 
low a more realistic polarization model. The work presented in 
this paper will be continued with such more detailed studies. 
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